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ADAPTIVE  CONTROLLERS* 
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and  Gunter  Stein 


Abstract 


aL 


There  have  been  many  algorithms  proposed  for  adaptive  control  which  will 
provide  globally  asymptotically  stable  controllers  if  some  stringent  conditions 
on  the  plant  are  met.  The  conditions  on  the  plant  cannot  be  met  in  practice  as 
all  plants  will  contain  high  frequency  unmodeled  dynamics  therefore,  blind  im¬ 
plementation  of  the  published  algorithms  can  lead  to  disastrous  results.  This 
paper  uses  a  linearization  analysis  of  a  non-linear  adaptive  controller  to  de¬ 
monstrate  analytically  the  following  design  guidelines  which  alleviate  some  of 
the  problems  associated  with  adaptive  control  in  the  presence  of  unmodeled  dy¬ 


namics  : 


1.  The  gain  in  the  estima* ion  mechanism  should  be  kept  small.  This  will 
make  the  estimation  loop  slow. 

2.  The  overall  effect  desired  of  the  controller  should  be  reasonable. 
Don't  try  to  make  the  algorithm  do  too  much. 

3.  The  system  should  be  sampled  slowly  to  alleviate  the  effects  of 
unmodeled  dynamics. 

The  points  made  are  further  demonstrated  bv  simulation  results. 
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1.  INTRODUCTION 

There  have  been  many  algorithms  proposed  for  adaptive  control  which  will 
provide  globally  asymptotically  stable  controllers  if  some  stringent  conditions 
on  the  plant  are  met.  The  conditions  on  the  plant  cannot  be  met  in  practice  as 
all  plants  will  contain  high  frequency  unmodeled  dynamics  therefore,  blind  im¬ 
plementation  of  the  published  algorithms  can  lead  to  disastrous  results  111. 
This  paper  uses  the  linearization  analysis  introduced  in  [21  of  a  non-linear 
adaptive  controller  to  demonstrate  analytically  the  following  design  guidelines 
which  alleviate  some  of  the  problems  associated  with  adaptive  control  in  the 
presence  of  unmodeled  dynamics: 

1.  The  gain  in  the  estimation  mechanism  should  be  kept  small.  This  will 
make  the  estimation  loop  slow. 

2.  The  overall  effect  desire  of  the  controller  should  be  reasonable. 

Don't  try  to  make  the  algorithm  do  too  much. 

3.  The  system  should  be  sampled  slowly  to  alleviate  the  effects  of 
unmodeled  dynamics. 

The  value  of  these  guidelines  is  demonstrated  by  examples  using  the  adap¬ 
tive  algorithms  of  Goodwin,  Ramadge  and  Cains  [3].  Further  analysis  using 
other  algorithms  is  performed  in  [4],  Note  that  while  the  guidelines  discussed 
in  this  paper  do  much  to  alleviate  some  problems  of  practical  adaptive  control, 
the  end  results  should  not  be  called  robust  because  the  problems  of  disturb¬ 
ances  in  adaptive  controllers  as  discussed  in  [4]  and  [5]  remain. 

2.  AN  ADAPTIVE  CONTROLLER  AND  THK  LINEARIZATION  ANALYSIS 

The  actual  plant  is  represented  by  the  equation: 

-d 

g-q  PB 

y(t)  -  -  |.(OI 


(!) 


-2- 


and  the  reference  model  is  given  by: 


yM(t) 


A(n) 


(r(t)l 


(2) 


M 

where  B,  A,  B*i,  A*  are  polynomials  in  q"l,  the  backward  shift  operator,  i.e., 
q“Mu(t)J  ■  u(t-l).  The  superscript  in  parenthesis ,  as  in  B^n\  gives  the 

n 

order  of  the  polynomial. 

The  system  ia  designed  assuming  dp*d  and  the  degrees  of  B  and  A  are  m  and 
n  respectively. 

Auxiliary  variables  are  simply  delayed  versions  of  the  input  and  output 
variables,  as  follows: 


wvi(t)  -  q“Mv(t)l 

i*0,l , . . . ,n-l 

(3) 

wui(t)  ■  [u(t)1 

i“0,l , . . . ,n-2 

(A) 

gwBi1"^ 

wr(t)  -  fr < t )  I 

A(n) 

M 

(3) 

The  scalar  control  input  to  the  plant  is: 

u(t)  ■  kT(t)w(t)  (A) 

with  k(t)  being  a  vector  of  time-varying  gains, 


The  specification  of  the  algorithm  is  completed  by  the  addition  of  the 
parameter  adjustment  mechanism: 


Jt(t)  ■  k p 


y  j£d(t)e(t) 

l-q-<>  lWr(t)«y(t) 
d 


00 


a.-,  nff***?  ■ 
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where 

e(t)  *  y(t>  -  y*(t)  ( 

»y(t)  ■  q“d[w(t)]  (1 

and  ygp  chosen  so  that  YRp  <2  (1 

The  controller  structure  is  displayed  in  Figure  1. 

RM^M 

Since  r(t)  is  prefiltered  hy  —  ■ —  the  goal  of  the  control  loop  is  to 

Si 

become  a  dead  beat  controller,  i.e.,  have  the  closed  loop  transfer  function 
•P —  ° — *-  ,  shown  in  the  box  in  Figure  1 ,  be  a  pure  delay  of  A  steps  with  a 


unity  gain.  From  Figure  1,  we  obtain 


u 

g*B*q~dP  k?Fpd  P* 


A(l-q-1Ku*(n_2>)-gpq  P  UK**"-1' 


If  the  order  of  the  plant  model  is  the  same  as  the  order  of  the  plant  and 
the  delay  is  modeled  correctly  there  is  enough  freedom  to  achieve 

“d— 

LI  *  -L  -  o“d  (13) 


and  the  controlled  plant  output  will  match  the  model  output. 


The  error  equation  is  given  below: 


e(t)  -  g-B*q  d  -  kT(t)w(t)  ♦  P  a>  -Afr(t)l  (14) 

The  error  system  is  displayed  in  Figure  2.  This  system  in  non-linear  as  w(t), 
w<j(t)  and  e(t)  all  depend  on  y(t). 

The  equations  of  DA2  are  analysed  by  linearising  the  system  about  a 
nominal  set  of  parameters.  The  signals  w(t)  from  eqns.  (3)-(5)  are  represented 


-d  /*.*-( d  „-<!) 


as  follows: 


w(t)  ■  w*(t)  *  6w(t) 


(15) 
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and  the  parameters,  k(t),  are  represented  as  follows: 

Jc(t)  -  k*  ♦  6k(t)  (1ft) 

Assune  that  both  the  parameters  and  the  outputs  are  close  to  their  desired 
values  so  that  6w(t)  and  £k(t)  are  small. 

Linearizing  the  system  about  w*(t),  k*  and  assuming  that  the  reference 
input,  r,  is  constant,  the  signal  w*  can  be  taken  as  a  constant  and  the  error 
system  of  Figure  2  becomes  the  linear  and  time-invariant  system  of  Figure  3. 


Define 


so  that  0  <  d*  <  1 . 


l+v*Tw* 


The  error  system  can  then  be  represented  as  in  Figure  4  with 


Gg(q"l) 


t*B*q“dP 


(l-q“d)k*A* 

r 


The  following  analysis  is  performed  by  first  choosing  a  desired  system 

*R*  ~dp 

* - -* -  ,  k*  about  which  to  analyze  the  system.  The  behavior  of  the  system 

A*  r 

near  the  desired  operating  point  is  established  hy  performing  a  root  locus  for 
the  error  system  of  Figure  4  using  d*  as  a  parameter. 

3.  ANALYSIS  OF  SYSTEM  WITH  NO  tfNMODF.l.FD  DYNAMICS 

Consider  first  the  case  when  the  system  is  modeled  properly.  As  was  men¬ 
tioned  in  Section  2  it  is  possible  to  achieve 

P  _  -A  / 1»\ 


k*  -  — 

r  *r> 
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Substituting  eqns.  (13)  and  (14)  into  eqn.  (1ft)  yields: 


The  d*-root  loci  for  the  error  system  of  Figure  4  with  Gg  as  in  eqn.  (20) 
are  given  in  Figure  3a  for  d*l  and  in  Figure  3h  for  d"3.  In  both  cases,  it  can 
be  seen  that  the  linearised  system  in  stable  if 

0  <  <  2  (21) 
Indeed,  this  in  exactly  the  condition  for  global  stability  given  in  [3]. 
Thus,  for  this  algorithm  with  no  unmodeled  dynamics,  the  local  stability  analy¬ 
sis  performed  in  this  section  gives  the  same  condition  as  the  global  stability 
proof  performed  in  [3], 

If  one  could  choose  y  so  that 

Y«l*gp  -  1  (22) 

the  error  system  itself  would  he  deadbeat,  i.e.,  the  error  would  become  sero 

after  dp  steps  where  dp  is  the  properly  modeled  pure  plant  delay. 

From  eqn.  (16),  one  can  see  that  d*  is  close  to  unity  for  a  large  range  of 

reference  inputs.  Therefore,  if  the  plant  gain  gp  is  known,  eqn.  (22)  can  he 

satisfied  for  a  large  range  of  reference  inputs  by  choosing. 

1 

T  “  ~  (23) 

ftp 

Indeed  the  choice  of  y  given  by  eqn.  (23)  is  implied  as  the  proper  way  to 
choose  the  adaptation  gain  in  (3].  Although  this  choice  of  gain  produces  a 
very  fast  and  well  behaved  error  system  when  no  unmodeled  dynamics  are  present, 
we  will  see  in  Section  5  that  the  gain  of  eqn.  (23)  will  be  far  too  large  in 
the  presence  of  unmodeled  dynamics.  Large  gains  in  the  presence  of  unmodeled 
dynamics  will  produce  unstable  adaptive  systems. 
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4.  AN  EXAMPLE  OF  ANALYSIS  IN  THE  PRESENCE  OF  1INMODELED  hYNAMICS 

When  there  ere  unmodeled  dynamics  present  in  the  plant,  the  desired  system 
can  no  longer  be  chosen  so  that  its  output  exactly  matches  the  model  output. 


The  desired  system  must  now  be  chosen  to  match  the  model  as  closely  as 
possible. 

The  example  used  in  this  section  will  consist  of  a  first  order  plant  with 


(24) 


(2S) 


In  order  to  obtain  a  discrete-time  system  which  is  equivalent  to  the  sys¬ 
tem  of  eqn.  (24)  the  standard  technique  of  discrete-time  control  system  analy¬ 
sis  called  hold  equivalence  was  used  (see  |6],  Section  3.4).  The  continuous¬ 
time  plant  of  eqn.  (24)  is  preceded  by  a  aero-order  hold  and  followed  by  an  im¬ 
pulse  sampler  which  is  synchronized  with  the  zero-order  hold.  The  resulting 
discrete-time  system  is  equivalent  to  the  original  continuous-time  system  in 
that  both  systems  will  produce  the  same  output  at  the  sampling  instants  if  the 
input  to  each  system  is  constant  from  one  sampling  instant  to  the  next. 

Although  anti-aliasing  filters  (see  (6)),  are  usually  included  in 
discrete-time  controller  designs,  such  filters  are  not  specifically  treated 
here.  Any  filter,  such  as  an  anti-aliasing  filter,  operating  upon  the  plant 
output  can  be  considered  as  part  of  the  plant.  Indeed,  since  the  presence  of 
an  anti-aliasing  filter  is  often  ignored  when  designing  the  adaptation  mechan¬ 
ism,  it  is  reasonable  to  consider  an  anti-aliasing  filter  as  part  of  the  unmod¬ 
eled  dynamics  of  the  plant. 
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For  our  initial  investigation  a  sampling  period  of  T"0.04  seconds  was 
used.  This  represents  fairly  fast  sampling,  since  it  is  approximately  ten 
times  as  fast  as  the  fastest  dynamics  in  the  plant.* 

The  discrete-time  description  of  the  plant  is: 

y(t)  .  (o.oo3M)(i.(i.m<r1Ki»2.7t:>,l-i)<ri _  <„, 

(1-0. 961q'1)< l-(0.547*i0.044)q-l)(l-(0. 547-10.044)0-! ) 

The  model  was  chosen  as  the  discrete-time  equivalent  of  eqn.  (25) 

yM(t)  -  — -2)q~1  fr(t)  1  (27) 

l-(.R8)q-1 

Note  the  presence  of  the  non-minimum  phase  aero  at  -2.763  in  the  discrete¬ 
time  version  of  the  plant  as  predicted  in  1 7 1 .  The  validity  of  the  linearized 
analysis  is  unaffected  by  this  non-minimum  phase  zero. 

5.  THE  IMPORTANCE  OF  THE  ADAPTATION  GAIN 

Figure  6  shows  a  k*-root  locus  which  determines  the  poles  of  the  nominal 

y 

controlled  plant  of  eqn.  (12)  for  the  plant  of  eqn.  (26). 

A  set  of  gains  which  will  produce  a  stable  nominal  controlled  system  is 
given  by: 

k*  -  -0.8;  k*  -  1.32  (28) 

y  r 

g*B*q"d 

The  position  of  the  poles  of  &  ”  ■  with  the  parameters  of  eqn.  (28) 

A* 

are  indicated  by  boxes  (O)  in  Figure  6. 

The  yd*-root  locus  of  the  error  system  poles  given  by  eqn.  (18),  using  the 
parameters  of  eqn.  (28),  is  shown  in  Figure  7.  From  Figure  7  it  can  be  seen 
that  the  error  system  is  unstable  for 

yd*  >  .15  (29) 


*The  importance  of  the  sampling  interval  in  determining  the 
stability  of  the  adaptive  control  system  in  the  presence  of 
unmodeled  dynamics  is  discussed  in  Section  7. 


;;  -iWgS**'  **■ 


Figure  8  shows  Che  results  of  s  simulation  generated  with  the  parameters 

initialized  at  the  values  given  in  eqn.  (28)  and  with 

Y-0.2;  r-0.1;  Yd*-0.04  (30) 

# 

The  system  is  well-behaved. 

Figure  9  shows  the  results  of  the  same  simulation  hut  with 

Y-0.2;  r-10.0;  y<**-.199  (31) 

As  expected  from  the  Yd*-root  locus  of  Figure  7  the  system  is 
oscillatory. 

If  yd*  is  too  large  then  there  exista  no  nominal  system  for  which  the 
error  system  is  stable.  This  is  shown  in  Figure  10  which  gives  the  k*-root 

y 

locus  of  the  linearized  error  system  of  Figure  4  with 

Yd*  -  .44  (32) 

For  this  numerical  example,  the  analysis  shows  that  if  yd*  is  close  to 
0.94,  instability  problems  will  ensue,  while  if  the  value  of  yd*  is  smaller 
that  of  eqn.  (32),  the  yd*-root  locus  of  Figure  10  will  pass  through  the  unit 
disk  indicating  that  a  set  of  parameters  for  which  the  linearized  system  is 
stable  exists. 

Thus  if  we  let  y*l  and  r  is  increased  yd*  will  approach  0.04  and  the 
system  will  go  unstable. 

Figure  11  shows  the  simulation  result  with 

Y*1.0;  r«1.9  (33) 

The  value  r*1.5  corresponds,  through  eqn.  (17),  to  d*".81.  In  this  case,  there 
are  sets  of  parameters  for  which  the  linearized  system  is  stable.  The 
parameters  of  the  simulation  converge  to  such  a  set  of  parameters. 

Figure  12  shows  the  results  of  a  simulation  with  r  increased  to 

r-3.l  (34) 

This  value  of  reference  input  corresponds,  through  eqn.  (17),  to  d**.95.  How 
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there  is  no  set  of  parameter  values  for  which  the  error  system  is  stable  and 
the  plant  output  blows  up. 

Thus,  the  first  design  guideline  has  been  shown: 

In  order  to  maintain  stability  in  the  presence  of  unmodeled 
dynamics,  it  is  necessary  that  the  adaptation  gain,  y,  of  the  system 
be  kept  small  and  the  adaptation  proceed  slowly. 

A  few  remarks  are  in  order  here: 

Remark  1:  The  adaptation  gain  of  around  1.0  which  is  the  limit  for  any  chance 
of  stability  is  an  order  of  magnitude  smaller  than  the  gain  which 
would  be  picked  using  the  guidelines  of  eqn.  (23)  and  ignoring  the 
high  frequency  unmodeled  dynamics. 

Remark  2:  While  arbitrarily  fast  adaptation  is  theoretical  possible  fft],  the 
above  example  demonstrates  that  the  concept  is  extremely  dubious 
given  more  practical  assumptions. 

Remark  3:  The  linearization  technique  shown  here  provides  good  guidelines  for 
the  limits  of  the  size  of  the  adaptive  gains  necessary  for 
stability.  Such  gains,  however,  may  not  be  sufficiently  small  for 
arbitrary  inputs  and  initial  parameter  estimates. 

7.  THE  IMPORTANCE  OF  THE  NOMINAL  CONTROL  LOOP 

In  order  to  study  the  ability  of  the  algorithm  to  match  the  reference 
model  in  the  presence  of  unmodeled  dynamics,  we  return  to  study  the  nominal 
controlled  system  of  eqn.  (12)  when  the  nominal  parameter  values  of  eqn.  (2ft) 
are  in  effect  along  with  the  plant  of  eqn.  (26)  and  the  model  of  eqn.  (27). 

The  nominal  controlled  system  with  the  parameters  of  eqn.  (2ft) 

g*B*q~d  _  (.0046)q-l(l*0.196q“l)(le2.763q-*)  (35) 

A*  (l-.82q“l)(l-.79q-1)(l-.A5q-1) 
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allows  the  nominal  closed-loop  controller  of  Figure  1, 


(t)  .  (. 0046) (.12)q“l( 1*0. 196q-l)( 1+2. 763q~*  [r(t)l 


(36) 


(l-.fi8q-1)(l-.R2q"1)(l-.79q-1)(l-.45Q“1) 
to  match  the  model's  d.c.  gain  hut  does  not  provide  an  overall  reasonable  match 


of  the  model.  In  fact,  no  values  for  the  paratmeters,  k*  and  k*,  will  allow 

Y  r 


even  approximate  model  matching. 

The  reason  that  this  algorithm  cannot  even  approximately  match  the  model 


in  the  presence  of  unmodeled  dynamics  is  that  the  nominal  control  loop, 
g*B*q~dP 


,  must  converge  to  a  dead-beat  controller  no  matter  what  the  model 


is.  Due  to  the  dead-heat  structure  of  the  algorithms  all  the  poles  of 


must  be  moved  near  the  origin  to  provide  for  good  model  matching. 


This  cannot  occur  for  this  example  as  Figure  6  shows.  The  prefilter  of  r(t)  by 
8m  ®m 


Am 


followed  by  a  dead-beat  controller  is  a  poorly  designed  nominal  control 


loop  in  the  presence  of  unmodeled  dynamics.  It  requires  much  larger  feedback 
gains  than  are  necessary  to  match  the  model.  Other  algorithms,  for  example, 

|9j  and  (in],  provide  structure  which,  if  chosen  correctly,  allow  approximate 
model  matching  in  the  presence  of  unmodeled  dynamics.  In  using  the  added 
flexibility  of  other  algorithms  the  lessor  of  this  example  should  remain. 

Design  the  nominal  control  loop  so  that  approximate  model  matching 
can  be  easily  attained  even  in  the  presence  of  unmodeled  dynamics. 
Remark  1 :  It  may  occur  that  the  dominant  dynamics  of  the  nominal  control  loop 
can  be  made  to  match  the  dominant  model  dynamics  but,  in  order  to  do 
so,  the  unmodeled  poles  are  moved  close  to  the  stability  boundary. 
This  situation  will  severely  limit  the  size  of  the  adaptation  gains 
for  which  such  a  nominal  system  is  stable  as  the  unmodeled  poles 
will  be  pushed  over  the  stability  boundary  in  the  d*-root  locus. 


'  *>■ 


~  *t-i  1 
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Remark  2:  The  linearization  technique  shown  here  ia  a  good  technique  to  use  in 
deciding  what  ia  a  reasonably  designed  control  loop  in  the  presence 
of  unmodeled  dynamics. 

8.  THE  IMPORTANCE  OF  THE  SAMPLING  INTERVAL 

There  ia  an  additional  design  parameter  in  discrete-time  systems  which  can 
increase  tolerance  to  unmodeled  dynamics.  That  parameter  is  the  sampling  in¬ 
terval  . 

Let  the  system  of  eqn.  (24)  he  sampled  at  a  rate  of  T“.4,  instead  of 
T*.04.  Such  a  sampling  rate  is  fairly  slow  in  that  it  represents  five  times 
the  speed  of  the  modeled  pole;  however,  sampling  of  the  unmodeled  dynamics 
occurs  at  barely  once  per  cycle.  The  equivalent  discrete-time  system  is  now 
given  by 

y(t,  .  _ (Q.«29)(l*.0319.-l)(l-..0W»q-l),-l _  |u(t))  (J7) 

(l-.67q-I)(I-(.n017*i.O01R)q'1 )(!-(. 001 7-i .OOIRJq'1 ) 

We  note  there  is  no  longer  a  non-minimum  phase  zero,  as  was  the  case  in 

eqn.  (26)  where  the  plant  was  sampled  faster.  Indeed,  the  poles  and  zeroes  of 

the  unmodeled  dynamics  are  very  close  to  each  other  so  that  their  effects 

almost  cancel. 

Figure  13  shows  the  k*-root  locus  of  the  nominal  controlled  plant  of  eqn. 

y 

(12)  when  the  open  loop  plant  is  described  by  eqn.  (37).  From  Figure  13,  we 
notice  that  all  the  nominal  control  system  poles  can  he  placed  close  to  the 
origin  and  that  the  nominal  closed-loop  controller  of  Figure  1  can  be  aMde  to 
match  the  model  fairly  closely. 

The  yd*-root  locus  of  Figure  14  then  shows  that  the  unmodeled  dynamics 
hardly  come  into  play  allowing  the  full  yd**2  of  eqn.  (21)  gain  with  retention 
of  stability. 


Figure  15  shows  Che  results  of  a  simulation  made  with  the  system  sampled 


at  T“.A  and 


r-10.0 


The  parameters  were  started  at  what  would  he  their  desired  values  if  no 
unmodeled  dynamics  were  present 

k*  -  -1.06;  k*  -  1 .58 

y  r 


(38) 


(39) 


The  system  behaves  as  if  there  were  no  unmodeled  dynamics  at  all.  The 
plant  and  the  model  output  in  Figure  15  coincide. 

Figure  16  shows  the  same  type  of  simulation  but  with  the  parameters  start¬ 
ed  out  at  zero.  The  system  adjusts  quickly  to  follow  the  model  again  as  if  no 
unmodeled  dynamics  were  present.  Thus,  the  third  design  guideline  has  been 
demonstrated: 

Sample  the  system  slowly  enough  to  remove  the  effects  of  unmodeled 
dynamics. 

Remark  1:  Anti-aliasing  filters  will  help  the  situation  if  they  are  considered 
part  of  the  modeled  dynamics  and  the  order  of  the  adaptive 
controller  is  increased  to  accommodate  the  filters.  Otherwise,  they 
simply  add  to  the  problems  of  unmodeled  dynamics. 

Remark  2:  Much  of  the  benefits  of  slow  sampling  depends  upon  the  movement  of 

the  zeroes  of  the  discrete-time  system  with  changing  sampling  rates. 
See  [7]  for  more  on  this  subject. 

Remark  3:  Again,  the  linearisation  technique  applied  here  provides  a  good 
final  check  to  see  if  the  sampling  is  slow  enough. 


9.  CONCLUSIONS 

In  this  paper  the  following  design  guidelines  for  adaptive  controllers 
have  been  demonstrated  analytically: 


1.  In  order  to  maintain  stability  in  the  presence  of  unmodeled  dynamics, 
it  is  necessary  that  the  adaptation  gain  of  the  system  he  kept  small 
and  that  the  adaptation  proceed  slowly. 

2.  Design  the  nominal  control  loop  so  that  approximate  model  matching  can 
be  easily  achieved  even  in  the  presence  of  unmodeled  dynamics. 

3.  Sample  the  system  slowly  enough  to  remove  the  effects  of  unmodeled 
dynamics . 

Although  the  demonstration  was  made  for  one  algorithm  via  example,  the 
guidelines  hold  far  more  generally.  In  f 4 ? ,  other  algorithm  are  used  in  analy¬ 
tical  demonstrations  of  the  guidelines. 

In  addition,  the  linearisation  analysis  introduced  in  [2]  and  used  in  this 
paper  is  shown  to  provide  a  method  of  analysing  adaptive  controllers  with  re¬ 
spect  to  how  well  they  follow  the  design  guidelines. 

Finally,  we  note  that  in  addition  to  the  problem.,  faced  by  adaptive  con¬ 
troller  which  are  alleviated  by  adhering  to  the  design  guidelines  addressed 
here,  there  are  other  problems  with  adaptive  controllers  (3)  which  need  to  be 
addressed  in  order  to  obtain  an  easily  usable  design  methodology. 
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